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ABSTRACT 


Hydrodynami cs of an axi -symmetri c combination bl owi 
steel making process has been analysed mathemati cal 1 y. Tht, 
reel rculatory flow field within the reactor vessel has beer 
represented interms of the turbulent Navi er -Stokes equations ir 
conjunction with the k ~£ tur bul ence model and the resultant set of 
partial differential equations together with the appropriate 
boundary conditions solved via a control volume based finite 
difference procedure embodied in the popular TEACH -T computer 
code. Prior to any evaluation of the model predi cti ons , sensi ti vi ty 
of computations to the choice of residual value, noda 
conf igurations and differencing schemes Ce. g. ,upwind, hybrid, 
exponential .power law etc. D have been assessed and results 
independent of these established. 

The bottom/submerged gas injection phenomena in a 
combined blowing reactor has been represented via the single phase 
variable density formulation technique, wherein the rising two 
phase gas-liquid mixture has been treated like a homogeneous fluid 
of reduced density. The gas volume fraction and hence the density 
of the two phase mixture have been estimated from a knowledge of 
the dimensions of the two phase reg: on and the latter *s average 
rise velocity. To this end, several calculation procedures 
C viz. , formulations based on slip vs no-slip, constant gas volume 
fraction vs variable gas volume fraction in the gas-liquid plume 
etc. D reported in the literature have been critically examined. 
Through extensive comparisons between theory and experiments it is 
demonstrated that the bottom injection phenomena can be 
mathematically adequately represented by assuming gas /bubble 
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slippage and considering a constant rise velocity of the two phase 
mi xtur e. 

In representing the impinging gas jet Ce. g. , the top gas 
injection system} the shearing force imparted by the outwardly 
flowing gas across the cavity surface Ce. g. * the gas*-, i qui d 
interface} has so far been assumed to be the only mechanism 
producing the bulk liquid recirculation. However , the present 
investigation appears to suggest that flows generated through such 
considerations alone are not sufficient enough to drive the flow 
recirculation observed experimentally. Consequently, as an 
additional mechanism to the shearing force* the force of 
impingement estimated from the increased potential energy of the 
melt Cvi2. , deduced from the corresponding increase in bath 
height} has been considered and the top blowing system 
represented accordingly via an approach based on these two 
mechanisms. Flow fields thus predicted have been evaluated 
against equivalent experimental measurements on water models and 
excellent agreement achieved. Finally* hydrodynamic models of top 
and bottom injection systems have been coupled and a steady state 
fluid flow model for combination blown steel making system 
proposed. Numerical predictions from the model have also been 
evaluated against experimental measurements and reasonably good 
agreement between theory and experiments 


demons tr at ed . 
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Coefficients of discretization equation 

r epr esent i ng the effect of convection and 
dif f usi on 

Summati on of al 1 nei ghbour poi nt 

coef f i ci ent s 


Constant defined in CI.4D: has a 

value of 7 

Dissipation rate constant : has a value 
of 0. 09 

Constant defined in Cl.©}: has a value 
of 0. 12 

Constant defined in Cl. 4}: has a value 
of 0.107 

Constant defined in CI.©}: has a value 
of 9. 1 

Constant defined in Cl I. 10}: has a value 
of 1 . 43 

Constant defined in C I I . 1 0} : has a value 
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Diameter of a lance nozzle, m 
Diameter of a gas bubble > m 
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Acceleration due to gravity, m s 
Hei ght of the 1 i qui d bath , m 
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base of the vessel , m 



Constant defined in CI.4D: has a value 
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Pressure within the liquid, Pa 
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Density of liquid* kg m 
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2 -3 
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Constant in the k ~£ model : has a value of 
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Chapter I 


INTRODUCTION 

1.1 Introduction To The Thesis 

The injection of reactive and non-reactive gases into 
melts contained in metallurgical furnaces, ladles etc. is a 
frequent phenomena encountered in numerous metal processing 
operations. In either case, the injected gas promotes fluid motion 
within the vessel and thereby exerts considerable influences on 
phenomena, such as, si ag- metal r eacti ons , thermal and particulate 
homogenisation, inclusion float out or re-entrainment etc.. 
Consequently, gas injection induced stirring in metallurgical 
reactors has since long been the subject of considerable interest. 
Moreover ,this has been recognized to be of vital importance in 
determining the efficiency of numerous associated processes 
Cviz. ,melting/dissolution, di sper si on^mi xi ng etc. } carried out in 
gas stirred reactors. 

Typically, two different modes of gas injection have 
been applied to produce the required stirring and these include : 
Ci} gas jet impinging on the surface of metallic melts 
and 

CiiDsubsurf ace injection of gas through a lance, tuyere or 

porous -pi ug. 

As an example to the first category, mention can be made of the 
top blown oxygen steelmaking systemC viz. , the LD process} while 
for the latter, depending on the intensity of the blow a wide 
range of metal processing operations Cviz. , from bottom blown 
steel making processes to ladle processing of steel} may be 
cited. It has been realised for some time that the intensity of 
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bath stirring in top blown reactors is relatively sluggish, as the 
transfer of momentum from the impinging gas stream to the bulk 
liquid phase is rather inefficient. On the contrary, bottom 
injection provides vigorous bath agitation and thus exacerbates 
the rates of various mass transfer controlled phenomena, 
characteristic of high temperature operations. Indeed, in 

combination blown steelmaking processes desired metallurgical 
results have been successfully obtained by judiciously combining 
the top C this supplies the necessary oxygen for metalloid 
oxidation reactions? and the bottom/submerged gas injection 
systems Csee Fig. I. ID. Since the primary objective of the bottom 
gas injection system in a combined blown reactor is to provide 
necessary stirring rather than supply oxygen C the latter is known 
to be detrimental, since the bottom/tuyere life decreases by 
injecting oxygen through the bottom, affecting considerably the 
turnover rates per converter? , consequentl y injection of inert 
gases like argon etc. have been typically practised. Furthermore, 
the "gas-liquid coupling" in submerged gas injection being more 
efficient Cthan an impinging gas jet systemD , r equires a small 
volume of gas to be injected in relation to the top blown gas 
volume C typically about 5%D to produce the required stirring 
ef f ect . 

Since the flow fields induced in combined blown reactors 
Csay, the Q-BQPD play paramount roles in controlling the 
associated processes, such as slag-metal reaction rates, 
melting/dissolution of scrap, refractory lining wear etc. 
consequently, adequate knowledge of bath hydrodynamics is critical 
to the overall process analysis. However, since high temperature, 
visual opacity and relatively large sizes of metallurgical 
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Refractory lined 
vessel wall 



Fig. 1.1. The combination blowing steel making 
system . 





reactors pose serious limitations to direct experimental 
observations, consequently, an alternative approach such as cold 
model studi es/mathemati cal modelling C or both in conjunction D 
has to be adopted to develop understanding of fluid flow and 
other associated phenomena. 

In the present work, hydrodynamics of a combination 
blown steel making process has been modelled mathematically and the 
model predictions assessed against experimental work reported in 
the literature on equivalent room temperature physical models 
Ce. g. , water model sD . Towards these, at the onset, a brief review 
of relevant literature has been presented and the scope of the 
present work discussed. In subsequent chapters, the development of 
the mathematical model and the numerical modelling techniques have 
been outlined in detail and finally, several computed results have 
been compared against experimental observations. Summary and 
conclusions as well as recommendations for further work appear 
respectively in the concluding chapters. Computer programs etc. 
have also been assembled in an Appendix. 
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I # 2 Review Of Relevant Literature 

During the last two decades or so numerous 

investigations have been carried out on the combination blown 
steel making process and phenomena such as, mixing, slag-metal 
reactions, mass transfer rates etc. have been extensively studied. 
However, fluid flow aspects of the combination blown process has 
remained relatively less well investigated. Furthermore, little or 
no effort has so far been made to address fluid dynamics of such 
systems from mathematical points of view. Consequently, barring 
only one instance, no studies on mathematical modelling of 
combination blown process have been reported in the literature. 

Similarly, mathematical mode] studies of the top gas injection 
system Cviz. , the LD configuration!) has also been rare and only 
one such study so far have been noted. In contrast, however, a 
large number of investigations on the mathematical modelling of 
submerged gas injection systems appear in the literature. Since 
combination blown process embodies both top and bottom injection 
systems, consequent! y , fluid model studies on these systems have 
also been included in the present review in addition to those on 
combination blown system. 

The bottom/submerged gas injection system in the 

combined blowing reactor involves the Injection of gas into the 

melt through one or a number of tuyeres located at the base of the 
vessel. Since the volume of gas injected through the bottom is 
typically small, the injection process consequentl y , leads to 
phenomena simila to those observed in ladle processing of steel. 
The injected gas in the immediate vicinity of the tuyeres rapidly 



devolves into plumes of spherical cap bubbi es * rise due to 
buoyancy , imparting energy to the bulk liquid and leaves the 
vessel through the free surface. Therefore, it is more appropriate 
to consider a buoyant plume rather than a jet model to describe 
the bottom injection phenomena. The interactions between the gas 
and the surrounding liquid induces flow recirculation in the 
vessel as has been mentioned already and this facilitates many 
metallurgical processes, such as, slag-metal reactions, 

dissolution, mixing etc.. 

There are two different computational procedures 
typically been adopted to mathematically describe the 
hydrodynamics of such two phase buoyancy driven systems and these 
i ncl ude: 

CaD models derived on basis of Eul er i an frame of 
reference, which are further classified into two groups viz. , 

C i 3 the two phase models 

CiiDthe quasi -single phase models 

and 

Cb3 models derived on the basis of Lagrangian frame of 

ref erence. 

In the two phase modelling, two sets of continuity and 

C 1 3 

momentum equations are solved, one for each phase .This approach 

% 

of modelling being relatively more complex from the computational 
view points, has therefore, received less attention in the 
literature. Momentum equations for individual phases are 
interlinked via the interface friction f actorCe. g. » slip3 and in 
addition, no assumption is made about the gas volume 

fraction C vi z . , voi dage3 distribution in the plume. The latter is 
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rather estimated through the solution of the gas phase continuity 
and momentum equations. 

The quasi —si ngl e phase modelling involves solution of 

only one set of continuity and momentum equations. Such an 

approach* however* requires a prior specification Ceither 

estimated or empi ri cal 3 of the plume geometry Cviz. .the gas- liquid 

zone) as well as the gas voidage. The buoyancy force per unit 

volume Cp^gcO estimated on the basis of the assigned gas voidage 

Cviz. * oO is added to the vertical direction momentum equation as 

the required driving force. Eulerian single phase modelling has 

received considerable attention in the literature because of its 

relative simplicity and ease of computational efforts. This 

techni que, f or example, has been used extensively by Debroy et 

al . , Szekely and coworkers , Sahai and Guthrie and 

C 7 0D 

Mazumdar and Guthrie > to describe turbulent recirculaiory flow 

fields in bubble stirred ladles. 

In the latter appr oach, Ce. g. *CbDD the second 

phased . e. * the gasD is treated essenti ally by a Lagrangian 

technique, in which a bubble trajectory equation is used in 

addition to the liquid phase continuity and momentum equati onsC the 

latter are same as those in the quasi- single phase modelling!). 

The bubble trajectory equation together with the momentum 

equations provide estimates of the two phase region, gas voidage 
C 9!) 

etc. 


Of the three approaches of modelling a submerged gas 
injection system mentioned above, the quasi -single phase modelling 
technique has been more popular. Consequently, the present 


1 i ter atur e revi ew 


the mathematical modelling of the 


on 
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hot t om ■-■'submer ged gas injection system has been restricted to the 

studies based on such an approach only. 

C 2D 

Debroy, Majumdar and Spa] dj ng "were the first to 

introduce the concept of a buoyancy driven plume in the 

metallurgical ladles. They considered the density deficit between 

the gas and the liquid phases to be the driving force Cviz. , the 

buoyancy] producing the flow recirculation and modelled this 

effect through adding the buoyancy force term to the axial 

C 2D 

direction momentum equation. The authors adopted the Eulerian 

single phase modelling approach Ci.e. .one set of continuity and 

momentum equationsD together with a bulk effective viscosity 

model C1C> ‘ > for describing the turbulent flow in the system. 

Symmetry conditions were applied to all the variables C si nee the 

configuration considered was axi -symmetric] . In order to estimate 

the density deficit Cor the void fraction] in the two phase 

region, zero slip condition Ci.e., the gas and the liquid move 

togetherD has been applied. Numerical predictions of flow 

variables were assessed against experimental measurements of 

CUD 

Szekely et al . and qualitative agreement between theory and 

experiments have been demonstrated. 

C3D 

El-kaddah and Szekely proposed a fluid flow based 

mathematical model for estimating sulphur partitioning between 
slag and the metal phases in a gas stirred ladle. The 

desulphurisation kinetics was described through prediction of 
turbulent flow fields via the solution of Navi er -St ok es equation 
in conjunction with the k —e model. The calculation procedure 


in addition involved prior specification of gas voidageCcO through 
the solution of an ordinary differential equation, viz. , 



. . . ci . 15 


dz 


P P p 
£2 CtanW u z 3 
z 


and incorporating in an expression, viz. , 


c*Cz3 = 


£Q -nr aCl-cO U 3 
g c s 


2n f C r U dr 

o J 


... C 1 . 25 


in which is the slip velocity. Fluid flow based predictions 

were compared with the observed rate of desul phur i sati on deduced 

from high temperature melts and reasonably good agreement between 

theoretical and experimental results has been reported. 

C43 

Grevei.Szekely and El -Kaddah subsequently reported an 

extensive fluid flow study in a water model of an argon stirred 

ladle. In this investigation, Eulerian single phase modelling 

together with an assumption of bubble slippage have been 

incorporated. The formulation of governing equations etc. were 

essentially identical to those reported earlier by El -Kaddah and 
C 33 

Szekel y . Excellent agreement between predicted and experimental 

velocity fields C deduced through Laser Doppler Vel oci meter D was 

achieved. However, predicted and experimental turbulence kinetic 

C 43 

energy were only in moderate agreement. The authors 

acknowledged that such discrepancies were perhaps associated with 

the incorrect selection of constants or with the inherent 

shortcomings of the k -s turbulence model for representing such 

C43 

recirculatory flow systems. Furthermore, the authors suggested 
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that the principal mechanism of momentum transfer in bubble driven 
systems is due to fluid convection, C rather than turbulent 
diffusion!) and thus predictions regarding the velocity fields were 
not expected to be very sensitive to the choice of a turbulence 
model . 

C 5 f3!) 

Sahai and Guthrie ’ through their extensive 

experimental and theoretical studies, postulated that bubbles 
forming at nozzles or orifices rapidly devolve into a plume of 
spherical cap bubbles and that this is true for many gas injection 

processes carried out in holding or transfer vessels. The rising 

bubbles impart energy to the surrounding liquid and produce flow 
recirculation within the vessel. Equating the total energy 

supplied by the bubbles to the turbulence energy dissipation 
losses in the bath, a simple macroscopic expression for estimating 
the average rise velocity of the bubble plume has been proposed 
according to : 


Q 1/3 L 1/4 


U = 4. 4 C D ... Cl . 3D 

P * 1/3 

In subsequent numerical prediction of flow etc. , they adopted the 

above relationship to estimate the gas voidage Cthey considered no 

slip condi tionD in the two phase region. Furthermore, fixed 

velocity at the axis of the vessel has been applied as the 

required boundary condition instead of the typical zero gradient 

conditions. Predicted values were compared with experimental 

measurements of their own as well as with those of Lehner et 
Cl 3D 

al . and reasonable agreement between theory and experiments 


has been demonstrated. 
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C 7 8Z> 

Mazumdar and Guthrie # proposed hydrodynamic models 

for various axisymmeirlc gas injection configurations found in the 

industrial ladle refining operations. The authors estimated gas 

void fraction via equation Cl. 3D. The hydrodynamic model in 

essence was based on the solution of turbulent Navi er -Stokes 

equation together with the k -c turbulence model. Due to some 

discrepancy between theory and predictions for the CAS 

C4D 

conf iguration the empirical constants in the k ~c model were 

changed by some adhoc values. Numerically predicted values Cflow 
and turbulence parameters) thus obtained were shown to be in 

excellent agreement with experimental results. 

C l 43 

Salcudean and coworkers proposed a three dimensional 

calculation procedure to describe non~axi symmetric gas injection 
conf i gurati ons through the solution of turbulent Navi er -Stokes 
equation, incorporating the k-e turbulence model. There were, 
however, large discrepancies between theoretical estimates and 
experimental measurements particularly in the immediate vicinity 
of the nozzle. Even though good agreement between predicted and 
experimental values were shown for elsewhere, the authors 
suggested that the discrepancy could be from numerical errors due 
to false diffusion phenomena Ce. g. , related to the use of the 
differencing scheme etc.D as well as due to the empirical 
constants employed in the k ~s turbel ence model. 



IZ 


In contrast to these, not much information is available 

in the literature on the mathematical modelling of top injection 

CIS 1 6D 

system. However » some studies * ~ report advantages of bottom 

injection over top injection processes but donot address these 

C 1 TO 

from f undamental /mathemati cal point of view. Asai and Szekely 

were first to propose a mathematical model for the top injection 

systems. The authors formulated the flow problem by the stream 

C 1 81 

f uncti on— vor ti ci ty based method together with the 

C 1 91 

Prandtl -Kol mogorov formula to calculate the turbulent 

C 171 

viscosity. Numerically predicted results gave qualitative but 

not quantitative agreement with experimental measurements reported 

earlier by Wakelin^^ 0 ^ on the CO^ -water sys terns. Even though the 

authors reported reasonable agreement between numerical 

predictions and experimental measurements. the assumptions 

incorporated in the model were paid very little or no attention 

Cviz. . the shape of the cavity was assumed to be known and the 

vor ti city considered spatially uniform, an adjustable parameter 

over the cavity surface etc.l. Furthermore, the mechanism leading 

to the transfer of momentum from the jet to the liquid bath was 

also not adequately considered and properly justified. 

C211 

Zhang et al . proposed a mathematical model for 

combined top and bottom blowing steel making process. The problem 

was formulated by statements of the turbulent transport equations 

for stream function and vorticity incorporating the k -e model. To 

describe the top injection process, the authors considered the 

momentum transfer from the impinging gas jet to the bulk liquid 

via shear stresses acting over the cavity surface CFig. I. ID 
C21 !) 

They further related the momentum flux Ci.e. , shear stress!) 
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across 'the cavity surface to the velocity gradient and the 
viscosity of the gas Csee also later D by an empirical relation : 


tC r 3 


4 C Q p 

a 9 9 r 2 

’ i/2 , 72 exp 1 ~ c ~!Tc ITT 3 3 


. . . CI . 43 


x o d 
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In which C , K and C are the empirical constants Cfor other 
a r 

symbols, refer to the nomenclature!). 

The associated bottom injection phenomena was modelled by the 
quasi -single phase modelling approach. However, the plume geometry 
was not assumed to be known C rather estimated!) and gas void 
fraction accordingly was deduced from a new empirical relation: 
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a exp 
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where a Is the voi d fraction at the axis of symmetry and is given 
m 

as; 

a CzD = C Q 2/3 t 3 1/3 CH-Z+h ...Cl. 63 

m t g O 

Cp L ~P g >g 

Numerically predicted flow patterns were compared with 
experimental measurements and qualitative agreement between the 


two has been demonstrated. 



1*3 Scope Of The Present Work 


The brief literature review presented in section 1.2 
clearly indicates that hydrodynamic modelling of combined top and 
bottom blowing steelmaking processes has practically received 
little or no attention. Since fluid dynamics play key roles in 
determining the overall process efficiency and that adequate 
informations on flow variables etc. in high temperature 
metallurgical reactors can be usefully derived via mathematical 
model studies, it is therefore, appropriate to direct some efforts 
in this direction so that some understanding of the h^rod^amics 
of a combination blown steel making process can be developed and 
used to draw useful inferences on associated phenomena such as, 
scrap dissolution, refractory wear, slag-metal reactions etc.. 

So far, as has been pointed out already, a mathematical 

model describing the turbulent recir cul atory flow in a combined 

C 21 D 

blown reactor vessel has been proposed by Zhang and coworkers 

The model predictions, however, have not been rigorously evaluated 

against equivalent experimental measurements. Experimental flow 

patterns have only been assessed against theoretical predictions 

and qualitative agreement between the two has been claimed. Since 

detailed flow measurement studies carried out on aqueous models of 

similar systems have been reported subsequently by investigators 

elsewhere, consequently the procedure adopted Cas well the 

C 21 D 

assumptions incorporated in the mathematical model D by Zhang 

et al . needs to be re-evaluated. Furthermore, Zhang and coworkers 
have also employed a new calculation procedure to model the bottom 
injection system that is much different from the procedures 
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adopted by the pr evious investigators' 12 83 . Since numerous 
calculation procedures can be applied to represent the* 
bot t om/s ubmer ged gas injection system* consequently, the procedure 
adopted by Zhang et al . as well as the suggested superiorly of 
thei r technique above those of others require critical 
examination. This follows since no rigorous comparison between 
theory and experiments have been made to justify the assertions 
made 

Finally* the investigation of Zhang and coworkers 

appears to suggest that the transfer of momentum from the 

impinging gas phase to the bulk liquid takes place only via the 

shear stress acting across the cavity surface. Consequently* the 

shearing force has been considered to be the only likely mechanism 

driving the flow field and accordingly* a bydr odynami c model based 

C 21 D 

on this concept has been proposed . However, no evidence has 
been put forward to justify such considerations. In view of the 
fact that numerous metallurgical systems are driven by potential 
energy Cviz. , filling ladles, bubble stirred ladles etc.}, the 
assumption of a " shear force driven flow " clearly warrants some 
additional considerations. It is with all these preceeding 
objectives in mind that the present research program has been 
initiated. To these ends relevant details have been summarised in 


the subsequent sections. 



Chapter II 


THEORETICAL CONSIDERATIONS 
II. 1 Assumptions In Modelling 

Combined top and bottom injection of gas into melts 

contained in industrial steel making units involve numerous processes 
that are inherently multidimensional in nature. The interactions among 
the various phases Cviz. , slag , metal and gas) present in the reactor 
vessel invariably lead to numerous complex phenomena, which include 
for example, waves and droplet formation at the free surface, thermal 
and particulate gradients in the melt, slag entrainment and so on. 
Consequently, to describe the gas injection induced flow in such 
multiphase and mul ti di mens! onal systems via an appropriate 

mathematical frame work, several idealisations have to be made and 

accordingly, numerous assumptions need to be invoked. 

In the present work, the fluid-flow model for the 

combination blown process has been derived incorporating the following 
assumptions,©, g. , 

C13 gas injection into a cylindrical vessel has been considered 
and the presence of upper slag phases in the system has been 
i gnor ed , 

Cii3 the wavy free surface at the top has been considered to be 
essentially flat, 

CiiiD the gas is injected through a nozzle located at the base of the 
vessel as well as a lance placed above the melt. The effect of 
multiple nozzles and tuyeres have therefore, been neglected. 
Furthermore, the top as well as bottom gas injection systems within 
the cylindrical reactor have been considered to be placed along the 
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symmetry axis. Consequently, cylindrical symmetry of the system has 
been assumed and hence, variation of the flow properties along the B 
direction has been ignored, 

Ci\0 the bubble plume formed as a result of the submerged gas 
injection has been idealised and its geometry and dimensions have been 
assumed to be known, i 

CvD the shape of the cavity at the free surface resulting from the 
impinging gas jet has been assumed to be parabolic and the dimensions 
of the latter known, I 

CviZ> the recirculating flow in the reactor has been considered to be J 

I i 

driven by the shearing forces imposed by the impinging gas jet and the j j 
buoyancy force exerted by the rising bubbles from the submerged gas 

II 

j njecti on system, 1 

I j 

CviiD the system is essentially isothermal 1 

I 

and finally, 

CviiiD the induced flow is turbulent in nature. j 

j 

A schematic of the idealised system used for the present j 

mathematical model representation is shown in Fig II. 1. 

II. 2 The Governing Equations 

Under the imposed assumption of cylindrical symmetry, the 

flow variables have been assumed to be invariant along the e 

co-ordinate and hence interms of the cylindrical polar co-ordinate 

system Cr ,e and , the statement of the problem may be conveniently 

formulated via a two dimensional elliptic: C©. g. , the? flow is 

C 2-8D 

recircul aiing) turbulent flow model Cviz. , a set of flow equations 

in conjunction with an appropriate turbulence model}. 
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II* 2. a The flow equations 

The flow equations in essence are the turbulent 
Navi er -Stokes equations Cvi*. ,an equation of continuity with two 
seper ate equations for z Cthe axial) and r Cthe radial} directional 
moment um/moti on) and are represented below : 

Equation of continuity 

du 1 aCrvZ) 

— + r= o . . . . CII . 1} 

#z r dr 


Equation of motion in axial direction 
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Equation of motion in radial direction 
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eff. 
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Cl 1 . 53 


dr 


F £ 's and F g in equations Cl 1. 33 and Cl 1. 53 are respectively the shear 
force and the buoyancy force per unit volume imposed by the top and 
bottom gas injection systems and are essentially the driving forces 
producing the flow recirculation in the combination blown reactor. 
These will be quantified interms of appropriate expressions in the 
subsequent sec ti ons . 

II- 2.2 The turbulence model 


The effective viscosity C p^3 , appearing in equations 

CII.23 through CII.53 follows from a turbulence model and in the present 

C*-$~03 

investigation, the popular and the most widely used ^ turbulence 

C 

model of Launder and Spalding C viz. , the k-e model 3 has been 

adopted. The governing transport equations for the turbulence kinetic 
energy, k, and its dissipation rate, e, can be represented interms of 
the cylindrical polar coordinates as: 

Equation of turbulence kinetic energy 
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in which, , the net source term can be represented as 
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Equation of dissipation rate- of turbulence kinetic energy 
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The auxilary relationships embodied in the above set of equations Cviz. 
equationsCII . 2D to Cl I. SOD are : 


P eff + 


CII. 11D 


P-j- ^*D ^ ^ ^ ^ 


CII. 12D 


and finally. 


p = a p G + Cl -cO p L CII. 1 3D 

The gas volume fraction, a, appearing in equation CII.13D is finite and 

non-zero only within the plume region Ci.e. ,<D < r < r D. Elsewhere in 

c 


CII. 



The 


f i vc 


the? domai n , Cr D r^D , a is zero and consequently, p—p 

I — < 

empi ri cal constants Cviz. , C . , C_ , G . g» and cr*. D embodied in 

± Cj p £ Jc 

turbulence model were assigned the values according 
recommendations of Launder and Spalding. These are summarised 


the k -£ 
to the 
i n Tabl e 
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II. 3 Modelling Of The Bottom/Submer ged gas Injection Process 

As mentioned in the preceeding sections, the bottom 

injection system in the combined blown reactor has been mathematical 1 y 

modelled via the quasi -single phase approach Cviz. , the variable 

density formulation! assuming that the dimensions of the rising two 

phase region Cviz. .the bubble plumeD is known C typically an 

empirically determined parameter!. The buoyancy force exerted by the 

rising bubbles within the two phase region acts vertically upwards 

Cviz. .along the axial di recti onD and consequently, F_, the buoyancy 

C 

force per unit vol umeC =p L goD has been added to the axial 

direction momentum equation Csee equation CII.33D.lt is instructive to 
note here that a, gas volume fraction is finite and non-zero only 
within the two phase region. It is through the prior prescription of 
the gas volume fraction, a. within the two phase region that the 
buoyancy force resulting from the submerged gas injection system has 
been modelled in the present work. 

The gas volume fraction within the plume can be estimated 
readily applying the principle of volume continuity and towards this, 
numerous alternative approaches exist Csee also latterD . Typically, 
the rise velocity of the gas-liquid mixture is estimated first and 
then from the knowledge of the plume dimensions, the gas volume 
fraction, a, is conveniently deduced via either. 
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2 n f C U r dr 
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. . . Cl I. 14D 


. . . . CII. 155 


if slip conditions are applied. 

In equations CII. 140 and CII. 153, r^ represents the radius of the two 
phase plume and the rise velocity of the gas-liquid mixture at any 
axial station, z. Further, U £ in equation CII. 155, is the slip velocity 
and is conventionally deduced from the knowledge of average bubble 
sizes in the system. 

The density field distribution Cviz. , p in equations CII. 25, 
CII. 35 etc. 5 in the flow domain has been derived on the basis of the 
continum approach Cviz. .equation CII. 1355 and the density field 
established a priori. This is illustrated Fig II. 2. 


II. 4 Modelling Of The Top Gas Injection System 

The configuration of the impinging gas system together with 
the resultant cavity formed on the melt free surface have already been 
illustrated schematically in Fig 1.1 and Fig. II. 1 respectively. To 

mathematically model the flow induced by such top gas injection 

Cl 75 _ . . C 21 5 

system, Asai and Szekely as well as Zhang and coworkers 

assumed that the geometry as well as the physical dimension of the 

cavity are known and further considered that shear stresses 



Boundary of the 
idealised two 
phase gas-liquid 


\ Gas ♦ 

\ liquid / 
l region / 


Three types of control volumes 

Control volume A (completely inside the two phase plume) 

p 0 - °Op L 

Control volume B (partly inside the two phase region and 
partly inside the bulk phase) 
p = V f p(1-0C) ♦ v b p L 

Control volume C (completely inside the bulk phase) 

P = P L 

Fig. II. 2. The continuum approach used to define the 
density distribution in the calculation domain 




transmitted across the cavity surface essentially translates into the 

bulk phase liquid recirculation. 

To quantify the shear stresses transmitted from the 

impinging gas to the bulk liquid interms of an appropriate expression, 

C SI 3 

the empirical formula suggested by Zhang and coworkers has been 

adopted in the present work and the variation of shear stress across 
the cavity surface represented according to: 
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in which 1 , d, Q represent the key injection parameters, while C , 

^ y a 

C , K are empirical constants Cviz. , see nomenclature!). On the basis 
r 

of equation CII.16!>, the axial and the radial components of the shear 

stress at various locations on the cavity surface have been estimated 

Cby resolving rCrl into two mutually perpendicular stress components!) 

and considering relevant area, shear layer thickness etc. , 

the corresponding specific volume shear forces Cviz. , F and F !> 

S , Z s , r 

estimated and added to the axial and radial direction momentum/motion 
equations as the representative body forces imposed by the top gas 
injection system Cviz., equations CII. 33 and CII. 533. 


To correlate F and F Cthe specific volume shear force 
s , r s , z 

components, see equations CII. 33 and CII. 533 with the shear stress 
tCt 3, the parabolic cavity surface has been idealised interms of 
several line segments as illustrated in Fig. II. 3. Furthermore, for 
each line segment, a char acteri sti c radial distance Cr3 has been 
chosen Cthis essentially corresponds to the distance of the geometric 
centre of the line segment from the axis of symmetr y3 and the shear 
stress estimated via equation CII. 1 63 . Assuming the shear stress to be 



Normal face area 
to r (r)„ , z' 


Normal face area 
to r (r) r , r' 


Axial component of shear force = z TirJ 2 
Radial component of shear force = r'r( r) r 


Fig. II. 3. Segmentised cavity contour and the procedure 
adopted to evaluate components of shear force / 
stress at any point on the cavity surface . 
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pr evalent over the entire line segment., the two mutually perpendicular 
stress components C viz. , along the z and the r axes respectively!! have 
been derived, multiplied by the cor respond! ng normal area Cviz. , see 
Fig. II -33 and added directly to the axial and the radial direction 
momentum equations as the required specific volume body forces C= 
F or F in equations CII.33 etc. 3. To this end an appropriate shear 
layer thickness has been assumed Csee also later3. 

II. 5 The Boundary Conditions 

The boundary conditions applied to the numerical solution of 
the set of partial differential equations Cviz. , equations CII.13 
through Cl I. 1033 are: 

Ci3 At the axis of the symmetry Cr=0 , O < z < L3 

v = O 

du dk Qe 

— O ; = O and = 0 

dr dr dr 

Cii3 At the free surface of the liquid Cz=L, O < r < R3 

u = O 

dv dk de 

s O ; = 0 and - O 

dz dz dz 

and at the side walls Cr = R, 0<z<L3 and bottom surface C z=0 , 0<r<R3 
u=0, v=0, k=0 and e-O 

Close to the vicinity of the rigid boundaries Cviz. , the walls and the 

bottom surface! where variations in the flow properties are steep, the 

momentum Cu and v3 and the scalar Ck and £3 fields have consequently 
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been modelled using the routine wall functions 


Ce. g. , logarithmic 



law of wall for the- parallel to wall velocity component etc. These 
have been illustrated via Fig. II. 4. 

II. 6 Numerical Solution Procedure 

Since F^, F^^ *nd r can be conveniently estimated from 

the procedures outlined already in sections II. 3 and II. 4, it is 
therefore readily apparent that the five principal differential 
equations Cviz. , equations CII.l}, CII.2}, CII.4}, CII.6Z> and CII.9}} 
together with their appropriate boundary conditions C viz. , secti on 
II. 5} can in principle be solved for the five dependent variables u. 
v, p, k and c. 

Thus to numerically solve the set of partial differential 

equations represented by equation 01.1} through 01.10} together with 

the associated boundary conditions and the auxilary relations, a 

C23} 

control volume based finite differnce procedure has been adopted 

in the present work. In this, the flow domain has been discretized 
into a number of non-overlapping control volumes, and the governing 
equations integrated around all such control volumes Cstaggered 
control volumes have been considered for the two velocity components} 
to result into a system of algebric equations, which for any general 
variable 0, may be represented according to: 

A p 0 P ■ Ve * V* + Vn + Vs + S c 011 • 173 

in which, A is the centre point coefficient Csee Fig. 11.5} while A , 
P & 

A w , A , A^ are the neighbour point C2 dimensional!) coef f i ci ents , 

W N d 

representing the combined influences of fluid convection and 

diffusion. Further, S represents the positive component of the source 

c 



- Free surface of liquid 

(z =1, 0<rSR] 

_ <5v _ <5k d£ 

u * 0 ' 6Z =0 ' 3z =0 Qnd 6z =0 


— Axis of symmetry 
(r = 0, 0^ z $ L) 

^jr = 0,vs0,^=0 and =0 


Inside wall — 
[r=R, 0$z$L) 
u =0, v = 0,k = 0 
and e s 0 


Wall functions applied 
in the immediate 
vicinity of wall 


- Bottom wall 
(z = 0, os: r^R) 

\u =0, v = 0,k = 0 and £ =0 


Fig. II. 4. Representation of the various boundary 

conditions applied to the numerical solution 
of the governing partial differential equations 
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term, Cviz. , of S^, etc. see equation Cl I. 3D for examples while 

the corresponding negative component has been embodied in the 
expression of , the centre point coefficient. The set of algebric 
equations thus derived have been solved via the Tri -diagonal Matrix 
Algorithm, incorporating a line by line solution technique. The 
nonlinearity and coupling involved in the algebric equations C these 
follow from the governing partial differential equations} have been 
taken care off by adopting to an iterative calculation procedure. 
Hence, a sequential rather than a simultaneous solution procedure has 
been embodied in the present calculation scheme. 

The pressure-velocity interlinkage in equations CII.1D 

through CII.5D have been solved by an implicit finite differnce 

C 24D 

procedure refered to as SIMPLE C Semi -Implicit Method for Pressure 

C 25} 

Linked Equations}. For analysis of the gas-liquid region, the GALA 
technique Cthe Gas And Liquid Analyser} was incorporated into the 
SIMPLE procedure. In this, the physical properties of the fluid 
mixture in a cell in the two phase region has been averaged on a 
volumetric basisCsee Fig. 11.2} This required the conventional mass 
continuity equation to be replaced by a volume continuity equation, 
Cviz. , equation C 1 1 . 1 }} such that volume of fluid entering a volume 
element equals the total volume of fluid leaving the control volume. 

Since an iterative calculation procedure has been adopted, a 
convergence criterion CNormalised Absolute Residual Source < O. 005} has 
been set on all the variables Csee also later}, and calculations were 
carried out till NAPS CNormalised Absolute Residual Source} value for 
u, v, k, £ and volume continuity all fell below this stipulated value. 



II. 7 The Computer Program 


The numerical solution procedure outlined in the preceeding 

section refers to a modified version of the popular TEACH -T computer 

C 83 

program, used extensively by Mazumdar to underscore the 

hydrodynamics of ladle gas injection operations. In the present 

C 83 

investigation, the modified version of the TEACH -T code has been 

applied incorporating several other additional modif i cations 
/improvements. These for example include, 

Ci3 adaptation to an axi symmetric submerged as well as to an impinging 

gas jet system through incorporati on of various force components Ce. g. t ! 

F , F and F into the computational procedure), 

B» s , z s , r 

l 

CiiD prescription of the density field in the flow domain via the GALA 3 

i 

J 

procedure, j 

I 

! 

CiiiZ) incorporation of different differencing schemes C upwind, 
exoponential etc.} for calculation of combined convection and 
diffusion fluxes Cviz. , the coefficients A^., etc. in equation II.17D 

and finally, 

C 

Civ3 a procedure C based on the cell porosity method 3 incorporated 
to exclude the cavity volume from the domain of calculation. 

The latter modifications will be discussed elaborately in the 
subsequent sections. The flow circuit of the computer program applied 
to the present investigation has been illustrared in Fig. II. 6. 












Chapter III 


RESULTS AMD DISCUSSION 

III.l Analysis Of Sources Of Error In The Computed Results 

To assess the validity of predicted results from the 

mathematical model outlined in the preceedi ng chapter, comparisons 
between theoretical estimates and experimental measurements are 

essential. However, prior to any comparisons it is important to 

ascertain that mini mum numerical inaccuracies are involved in such 
calculated results. For a given set of equations and boundary 

conditions, the extent of assosiated error in the computed results is 
expected to depend upon several numerical parameters. These include, 

CaD the maximum allowable residual source Ce. g. ,the convergence 
criterion) , 

CbD the distribution of grid systems in the flow domain 
and 

CcD the differencing schemes applied to estimate the convection and 
diffusion fluxes embodied in the discretization equations Ce. g. , the 
coefficients A^, A^, etc.}. 

It is important to recognise here that these numerical parameters 
simultaneously affect the accuracy of results as well as the 
assosiated computation times, CPU. Consequently, an optimum balance 
between accuracy and CPU time has to be considered for any effective 
error analysis. 

III. 1.1. The influence of residual source sum values on computed result 

Since an iterative calculation procedure has been adopted to 
numerically solve the governing partial differential equations, a 
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convergence criterion needs to be specified. To this end, a residual 
source, R^. defined as 

R L A nb 0 nb + S c " A p 0 p l --CIII.1D 

has been estimated at each nodal point and added up over the entire 

flow domain for each individual variable involved Cviz. , u,v,K,£ 

etc. 3. Obviously, when the discretization or the finite difference 

equations are completely satisfied, F? L will be zero for all equations. 

However, since this requires an extremely large computation time to be 

achieved, a suitable convergence criterion is required to be specified 

so that the value °f ^ ^ J indicates summation over the two 

dimensional flow domain) be less than a certain small number. However, 

smaller is this allowable residual source value, the larger would be 

the computation time. Hence, it is desirable to find out an optimum 

allowable residual source value for the present problem. 

To assess the sensitivity of computed results to the value 

of residual source, R^, calculations were carried out for a submerged 

qas injection configuration C F_ is finite, F and F are 

ts s , z s , r 

respectively zero in equations CII.33 and CII.533 corresponding to the 

C 273 -5 

experimental conditions of Mazumdar et al . Ce. g. , Q=1 . 667x10 

3 

m Xs3 . These are shown in Table III.l for ready reference. A 

relatively large number of grid system Cviz. ,15x103 have been applied 

to encompass the flow domain CL=0. 21m and R=0. 15m) providing an 

effective grid spacing of approximately 20mm in either direction. 

C 233 

Further, the hybrid differencing scheme has been applied to 

estimate the convection and the diffusion fluxes Cviz. , the 

coefficients A , A., etc. 3. 

p N 
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Table III.l Experimental conditions of different investigators used in 
the computation of results presented in Figs. III.l through III. 5 



Experi mental 

conditions of 



Mazumdar and 

. C27D 
Guthrie 

Cast! 1 1 ej os and 

Br i mac ombe^ 

Lehner and co~ 
Cl 

workers ^ 

Bath height „m 

O. 21 

O. 4 

1 . 08 

Radius of the 
vessel * m 

0. 15 

0. 25 

0. 76 

Diameter of the 
nozzle , m 

1. 5x1 0~ 3 

6. 35x1 CT 3 


Gas flow rale, 

m /s 

Cl . 67x1 0" 5 , 

3. 33x1 0~ 5 
and 5. 0x1 0~S 

3. 6395x1 0~ 4 

1 . 667x1 0~ 3 1 

1 

1 

| 

i 

j 
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The predicted radial distribution of vertical velocity 

component at two different depths in the water model has been shown in 

Fig. III-l for three different values of ^ ^ Cviz. ,0.001, 0.005 and 

0. 025} . This appears to suggest that a reasonable value for the sum of 

the residual source is approximately 5xlO _3 and further lowering of 

this value does not lead to any improved accuracy in the predicted 

results. However, if a sum of residual source value of approximately 

0.025 is applied in the numerical calculation procedure, significant 

accuracy in predicted results appears to be lost. Consequently, for 

all subsequent calculations reported in this thesis, an allowable sum 

of residual value equivalent to 5x1 0~ 3 has been applied. Towards this, 

it is important to note here that in numerous equivalent previous 

studies C see for example reference 45, the selection of an appropriate 

convergence criterion has neither been adequately considered nor 

explicitly specified. In addition, usage of widely different 

C 7 21 5 

convergence criteria has also been noted ’ . This is particularly 

critical when comparisons between theory and experiments are to be 
made, since computed results depend considerably on such criterion. 

III. 1.2 The influence of nodal configurations on computed results 

The derivation of the discretization or the finite 
difference equations from the governing differential equations require 
numerous assumptions to be invoked Cviz. ,the mode of variation of the 
dependent variables between two nodal points etc. 5, and consequently, 
predicted results, as one would normally anticipate, would be 
sensitive to the nodal configurations applied to the numerical solution 
scheme. However, it is well known that in the limit of a very large 
number of grid network, these assumptions and their resultant 
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Dimensionless radial distance, 

Fig. III.1. Sensitivity of computations to the choice of 
residual (NARS) values illustrating the 
radial variation of the vertical component 
of flow for_ three typical NARS values 
(Q = 1.67x10' 5 m 3 /s, R=0.15m, L = 0.21m, 

15x10 grid system ). 
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influence on calculated results tend to be Insignificant and 
therefore, the numerical predictions appear to approach to the exact 
solution Cin situations where such solutions can be derived) . 
Consequently, it is essential to establish grid independent solutions 
prior to any comparison between theory and experiments. Towards 
this, clearly an optimum nodal configuration for each individual 
problem has to be known C Indeed, a more meaningful parameter would be 
the grid spacing rather than the numb&r of grids, since the former 
would allow one to deduce the required number of grids that would lead 
to grid independent results for any domain sizeD. 

To test these, flows in a submerged gas injection system in 

a cylindrical ladle CL=0. Sim and R=0. 1 5m) has been modelled 

mathematically by applying three different grid conf i gurati ons viz., 

10x7, 15x10 and 17x12 grid systems respectively. While, the hybrid 

C 233 

dif ferrencing scheme has been applied , an allowable sum of 

residual source equal to 0.005 has been considered in the calculation 
procedure. A comparison among numerically predicted radial 

distribution of vertical velocity component at two different depths in 
the water model is shown in Fig. III. 2. This indicates that both 15x10 
and 17x12 grid systems produce results that are practically 
equivalent. Consequently, predictions via a 15x10 grid system can be 
considered to be essentially grid independent. On the basis of this, 
grid spacing of approximately 2Gmm have been applied in eibher 
directions in all subsequent computations. It is to be noted here that 
the 10x7 grid system leads to results considerably different to those 
deduced from the two other nodal conf i gurati ons. 



Vertical velocity component , (m/s) x 10 



Fig. III. 2. Sensitivity of computations to the choice of 
nodal configurations illustrating the radial 
variation of the vertical component of flow 
for three different grid systems (Q =1.67x10 
rrr/s, R = 0.15 m, L = 0.21 m , NARS = 5x10 ) 



42 


III. 1.3 The influence of differencing schemes on computed results 

As mentioned in the previous chapter, the discretization 
equations embody coefficients CA p> A^, A s etc. > that represent the 
combined influences of fluid convection and diffusion processes. To 
estimate these coefficients Ce.g. , the convection and the diffusion 
fluxes) from the relevant velocity and diffusion coefficient values, 
several schemes can in principle be applied and these, for example, 
include : 

C1D the central differencing scheme 
C23 the upwind differencing scheme 
C33 the hybrid differencing scheme 
C43 the exponential differencing scheme 
and 

C 53 the power law differencing scheme 

It is to be mentioned here that except for the central differencing 

scheme, all others are known to produce physically realistic solution 

C23D 

for both high and low Reynolds number flows . Consequently, the 

central differencing scheme has not been considered in the present 

analysis because of such obvious limitations. 

Since these differencing schemes incorporate various 

assumptions Cfor example, the upwind scheme ignores diffusion effects 

C 23D 

for all values of Peclet number etc. D , consequently, numerically 

predicted results are likely to be a function of the differencing 

sch&m&s employed. Thus to test the sensitivity of computed results to 

the choice of the differencing schemes, flow calculations were 

performed via all the four schemes corresponding to the experimental 

C 27D 

conditions of Mazumdar et al. Cviz. . those used in the computation 
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of Figs. 111*1 and III. 2 Csee Table III. I}} and the predicted radial 
distribution of vertical velocity component at two different depths in 
the model has been compared in Fig. III. 3. This clearly indicates that 
the differencing schemes employed have practically negligible 
influence on the overall accuracy of predicted results. It is to be 
mentioned here that a 15x10 grid system together with an allowable 
residual source sum of 0.005 have been applied to derive the results 
presented in Fig. III. 3. 

Figure III. 3 further illustrates that diffusive momentum 

transfer in such systems is relatively less important * since results 

of the upwind and the rest of the differencing schemes are very 

similar. Such computed results therefore, appear to indicate that 

flows in metal 1 ur gi cal ladles are largely dominated by the inertial 

forces and that viscous effects Cboth molecular and turbulent} play 

insignificant roles in affecting the overall flow pattern induced. 

This is clearly substantiated from the results presented in Table 

III. 2, which indicate that the predicted mean speed of liquid 

recirculation via all the four differencing schemes are practically 

identical C maxi mum variation ± 15*0 and eqi valent to those observed 

C 

experimentally by Mazumdar et al . 






Fig. III.3. Sensitivity of numerically predicted vertical 

velocity component in a water model of a gas 
stirred ladle 77 illustrating the influence of diffe- 
rencing scheme employed on predicted results. 




Table III. 2 
via various 
measur emenis 


Theoretically predicted mean speed of liquid recircul ation 

differencing schemes and thei r comparison with experimental 

C 27D 

of Mazumdar et al . 


Differencing scheme 


Mean speed of liquid r eci r cul ati onC mm/sD 


Numerical predictions Experimental values 


Upwi nd 

37. 

13 


Hybrid 

40. 

07 


Power law 

40. 

07 

34. 0 


Exponent! al 


45. 3E 
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III. 2 Evaluation Of Calculation Procedures For Modelling Of The Bottom 
Gas Injection Process 

In the present study, as has been mentioned already in the 
preceeding chapters, the Eulerian quasi -single phase approach has been 
adopted to model the bottom gas injection process. Because of the 
relative ease with which such computational procedure can be applied, 
the quasi -single phase modelling has been widely popular and 
therefore, several versions of the procedure are available in the 
literature. In essence, the quasi -single phase modelling procedure is 
applied via a prior specification of the geometry as well as the gas 
volume fraction, a, in the plume. Since the gas volume fraction 
specified within the plume region plays key role in determining the 
overall flow field in such buoyancy driven systems, therefore, the 
various procedures available for gas volume fraction calculation are 
required to be carefully examined. To this end, any comparison among 
these calculation procedures has so far not been attempted, and 
consequently it is necessary to assess the adequacy of these by 
evaluating their relative performances, so that the most accurate 
approach can be identified and adopted to model the submerged gas 
injection system in a combination blown reactor. 

III. 2.1 Prediction of flow parameters and average speed of bath 
recirculation via different procedures and their comparisons with 
experimental measurements 

To estimate the gas volume fraction, a, in the rising two 
phase plume, different group of investigators have adopted widely 
varying concepts and based on these, four distinct types of 



computational procedures have so far emerged. These are discussed 
below in detail. 

r 3 4 ) 

III. 2* l.a The procedure adopted by Szekely and coworkers * 


Szekely and 
small volume element 
di f f er enti al equati on 
velocity viz. , 


C3* 4D 

coworkers through momentum balance over a 

in the plume region proposed an ordinary 
for estimation of the aver age plume rise 


dU 

z 


dz 


U 


z 


z 


71 Q g 9 

...cm. 13 

_ . _ 2 2 2 
2 Ctan X3 u z 
z 


in which, U z is the average rise velocity at any axial station, z. 

The boundary conditions considered applicable to equationCIII . 13 is 


The average rise velocity has been further related to the gas voidage 
a according to: 


Q - nr 2 aCl -cO U 
g c s 

aCz3 = . . .CIII. 2D 

r 

2tt f C U r dr 

o J z 


in which, is the slip velocity and is considered to be synonymus to 

the rise velocity of a characteristic single bubble in the qui scent 
liquid C5:> and is deduced from a knowledge of the average bubble 
diameter in the system viz. , 

gd B 1/2 

U aU B ■ 1.09C ~ -- 3 CIII. 33 

Further, U Q in the boundary condition is the superficial gas velocity 


and is deduced from : 
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n r 


CIII . 43 


The void fractions, a , 'thus obtained are clearly cross-sectional 
average void Tractions. It is apparent that once the plume geometry is 
known and the void fraction distribution within the plume established 
via equations CIII. 13 and CIII . 23 , the equations of motion can be 
readily solved to yield the distribution of flow variables in the 
entire calculation domain. 


C S 03 

III. 2. l.b The procedure adopted by Sahai and Guthrie * 

Sahai and Guthrie used an alternative concept to estimate 

the volume fraction distribution in the rising plume. Through detailed 

C 5 03 

hydrodynamic considerations, the authors * suggested a simple, yet 
effective formula for estimation of rise velocity in the two phase 
plume according to : 

A. S3 1/4 
Q L 

U = 4.4 CIII. S3 

P 1 


and related the plume rise velocity to the gas volume fraction via 
the volume continuity principle : 


Q 


U 


av 


2 . 
fi r L 
av 


CIII. 63 


in which, r refers to average radius of the plume, which is equal to 
av 

C 1 /*/33 r e y e Ce. g. , the corresponding radius at the mid bath depth3 . In 
equation CIII. 63 no consi derati on of bubble slippage has been given 
and further, the gas volume fraction has been considered to be 



spatially uniform over the entire volume of the two phase plume. 


III. 2. l.c The procedure adopted by Mazumdar and Guthrie 


C 281 


C p0D 

Mazumdar and Guthrie modified the procedure adopted by 

C 0 , 6> 

Sahai and Guthrie ’ ' to estimate the gas volume fraction in -the 


plume. These authors 


C £8} 


adopted the macroscopic plume model of Sahai 


and Guthrie Cviz. , equation CIII.4DD, and deduced -the gas volume 
fraction in the plume from the following relationship : 


aCzD = 


Q -nr cxCl-cO U 
g c s 


2 n U f C r dr 

D J 


.... cm . 5} 


Clearly, the value of is constant throughout the two phase plume 

and therefore a, the cross-sectional average void fraction is a 

function of z only. Equation CIII.03 is different from equation 

CIII.2} since instead of U , U , a constant plume rise velocity has 

z p 

been used. The authors justified the use of a constant rise 

velocity rather than a variable one Cviz. , U 3 in equation CIII.SD, 
J z 

C 1 3D 

since experimental plume rise velocity deduced by Lehner et al . 

were found to be independent of the axial distance. The procedures 

C S 0 28D 

adopted by Guthrie and coworkers * * like those of Szekely et 

al.^^*^ also require a prior specif ication of the plume geometry. 


r 04 ) 

III. 2. l.d Zhang and coworkers formulation for estimation of gas 


void fraction in the plume 


C PI 

Zhang et al . proposed a new approach to estimate the gas 

volume fraction within the plume. For the three previous approaches 
outlined earlier » the geometry of the plume has to be specified cl 
priori. In contrast, however » Zhang et al . * s formulation does not 



require prior specification of the two phase domain and instead, the 
gas volume fraction within the flow domain is deduced from the 
following empirical relationship: 


nCr.z) = a exp [ - — — 

m _2 


1 ...C1II.7D 


C CH-z+h_J 
9 O 


in which, ct represents the void fraction at the axis of symmetry and 
is estimated from : 


2/3 „ HL ,1/3 -3/5 


a Cz) = C. Q I =- 

m t Cp - p r?' 3 

y 


3 CH-z+h o :> . . C 1 1 1 . 83 


In equations Cl II. 73 and CI1I. 83, C , represents two empirical 

dimensionless constants which were assigned to values 0.12 znd 9.1 
respect! vel y. 

To assess the relative accuracy of predictions via the four 
different approaches outlined above as well as their adequacy with 
reference to the mathematical modelling of submerged gas injection 
systems, flow calculations were accordingly carried out corresponding 
to the experimental conditions reported earlier by Mazumdar et al . 
Csee Table III. 13. 

Towards this, to solve equation CIII.13 and hence to deduce 

the gas volume fraction in the plume a computer program incorporating 

C293 

a fourth order Runga-Kutta method has also been developed. This is 

presented in the Appendix. The plume dimensions were estimated from the 
previously published experimental data and further a bubble diameter 
of 1 Omm has been considered to deduce the slip velocity. Incorporating 
the turbulence model with each of the previously mentioned gas 

volume fraction calculation procedures, the distribution of the axial 
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CuD and the radial Cv) velocity components in the flow domain were 
deduced through the solution of relevant momentum balance equations 
and from these, the average or mean speed of liquid recirculation, 0, 
estimated as : 

J J* | Cu^+v^}*^! 2nr dr dz 
0 0 

u = — — . ..cm. 9} 

L R 

Bnr dr dz 



The above integration has been carried out numerically adopting the 
Eulerian integration procedureC vi z. , assuming to prevail 

over a control vol uim?) . 

Numerically predicted mean speed of liquid recirculation 

thus obtained have been directly compared with the corresponding 

C27D 

experimental measurements of Mazumdar et al . derived for three 

different gas flow rates in a water model C L=0. 21 m and R=0. ISnD of a 
gas stirred ladle system. It is to be mentioned here that the kinetic 
energy of the incoming gas has been ignored for the flow calculation 
as a first approximation, since this is significantly lower C typically 
about lO % or less) in comparison to the potential energy afforded by 
the gas rising through the liquid. 


The comparison between theory and experiments illustrated in 
Table III. 3. clearly indicates that the hydrodynami cs of the submerged 
gas injection system is adequately represented via the procedure 
suggested by Mazumdar and Guthrie . It is further apparent that 
Sahai and Guthrie’s prescription of gas voidage leads to a higher rate 
of bulk fluid recirculation. This can be anticipated since in their 
calculation procedure^’ 65 bubble slippage has not been considered. 

CENTRA ! ."~'RARY 

i' . ■ * 

"« 1 09084 . 
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III. 3 Experimentally measured 


speed of liquid 


, . . C 373 , , , . 

recirculation and their comparison with those estimated 
theoretically to illustrate the sensitivity of prediction to the 
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predicted mean 
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speed of bath 
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55. 66 


71 . 95 
15. 29 


53. 42 



With reference to the results shown in Table III. 2 it is also 
important to recognize here that Zhang et al ’ s formulation, in 
contrast to the three other equivalent computational procedures, 
provides significantly lower estimates of mean speed of liquid 
recirculation. 

III. 2. 2 Estimates of gas volume fractions within the plume and their 

comparison with experimental measurements of Castillejos and 

„ . . C3<X> 

Brimacombe 

Extensive measurements of gas volume fractions within the 

rising two phase plume in water models of gas stirred ladles have been 

C 

recently reported by Castillejos and Brimacombe ” . To assess the 

adequacy of the theoretically predicted voidage values deduced via the 
four different procedures outlined in the preceeding section, the gas 
volume fraction distribution within the rising plume were estimated 
corresponding to the experimental measurements of Castillejos and 
Brimacombe Csee Table III. 2D. 

The experimental centreline gas voidages have been compared 

directly in Fig III. 4 against the corresponding estimates derived via 

each of the four calculation procedures assuming a bubble slip 

velocity of 0. 22m/s Cthis corresponds to a bubble diameter of iOmmD. 

There, it is clearly seen that the procedure suggested by Mazumdar and 

Guthrie C28:> provides the most reasonable estimates of the gas volume 

fraction distribution in the plume. Similarly, in Fig III. 5 the 

cross-secti onal average void fractions deduced via. i^he four 

calculation procedures have been compared with equivalent experi mental 

observations of Castillejos and Brimacombe . Once aqain adequacy of 

C 282 

the procedure suggested by Mazumdar and Guthrie is readily 



Centre line gas vodage 


• Experimental measurements of 
Castiliejos and Brimacombe^ 

Following procedure of : 

— Mazun.dar and Guthrie^® 

— Sahai and Guthrie^ 

— Szekely et al.^ 



Fig. III.4. Theoretically estimated centre line gas voidages 
in a water model ladle and their comparison 
with equivalent experimental measurements 



Cross sectional average gas volume fraction 
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apparent. 


Referring to Figs. III. 4 and III. 5. one also finds that the 

gas volume fraction deduced via the procedure of Sahai and 

Guthrie is somewhat over estimated along almost 75% of the plume 

height. Similarly, the void fraction deduced via Szekely et al. *s 
C3,43 , 

procedure is somewhat smaller than those experimentally observed. 

In contrast to all these, Zhang et al . *s formula*" 21 "* provides 

estimates of gas voidages that are significantly lower than the 

experimental values. More important, however, is their*" 21 "* estimates 

of centreline gas voidages that seem to increase with the axial 

distance, z. which clearly contradicts the experimental findings. 

On the basis of the preceeding discussion, one would 

therefore expect that the hydr odynami cs of a submerged gas injection 

system can be simulated adequately via the calculation procedure 

C *2©D 

suggested by Mazumdar and Guthrie . While one would also expect 
that Sahai and Guthrie’s prescription of a would provide higher rates 
of bulk recirculation, equivalent procedures of Szekely and 
cowor kers*" 3 ’ ^ as well as that of Zhang and coworkers^" 5 would 
produce smaller rates of bulk recirculation in the reactor vessel. 
Such a trend in predicted results is readily apparent in Table III. 2. 
It is instructive to note here that in all the calculation procedures 
adopted for the estimation of a, barring the procedure of Zhang and 
coworkers , radial variation of gas voidages in the plume has been 
ignored*" 3,4 ’ 5 ' 6 ’ 28 "* and hence the cross-sectional average void 
fraction at any axial distance can be visualised to be identical to 
the gas voidage at the axis of symmetry Cviz. , the centre lineD. The 
results presented so far also appear to indicate that such radial 
variations can be ignored as a first approximation, so far as the 
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modelling of the bulk phase hydrodynamics is concerned. 


III. 2. 3 Prediction of rise velocity within the plume and their 
comparison with experimental measurements of Lehner et 

far, adequacy of the calculation procedure of Mazumdar 
„ . . . C 285 . 

and Guthrie has been demonstrated with reference to the estimation 

of gas voidage distribution in the plume as well as the mean speed of 

bath recirculation. To assess the appropriateness of the calculation 
C 2SD 

procedure further, particularly with reference to the prediction 

of rise velocity within the plume region, computations were also 

carried out corresponding to the experimental conditions of Lehner et 
C 1 3D 

al. C these are summarised in Table III. 15 and the predicted rise 

velocity at various axial stations within the bath compared directly 

C 1 3D 

in Fig. JII.fi with the corresponding experimental measurements 

Reasonable agreement between theory and experiments is once again 

readily apparent. With refernce to the results presented in Fig. 

III. 6, it is to be mentioned here that a slip velocity CU s 5 equivalent 

to O. 4 m/s has been assumed. This corresponds to a bubble diameter of 

approximately 3<5mm, which is typical of such water model systems in 

the regime of the gas flow rates considered. 

On the basis of the proceeding discussion, therefore, the 

C 28D 

calculation procedure adopted by Mazumdar and Guthrie has been 

applied to represent the bottom gas injection process in the 
combination blown system. 



Dimensionless radial distance, 



• Experimental measurements of 
Lehner et al 

Present calculation applying the 

procedure of Mazumdar and Guthrie^ 

Fig. III. 6. Theoretically predicted vertical velocity 
component within the rising plume and 
their comparison with experimental 
measurements of Lehner et al. 




HI* 3 Modelling Of The Top Gas Injection System And Comparison with 
Experimental Measurements 

In the preceeding sections, modelling of the bottom gas 
injection process has been considered in detail and an adequate 
mathematical representation for the gas injection induced flow 
developed. In this section, modelling approaches for the impinging gas 
jet system currently available in the literature will be discussed and 
evaluated rigorously. Once the appropriate procedure has been 
identified, the desired mathematical model for the combination blown 
system can then be derived by considering these models for the bottom 
as well as the top gas injection systems respectively. 

III. 3.1 The shearing forces acting across the cavity surface and its 
role In driving the bulk phase flow In the top gas Injection system 

C 1 7 21 3 

Earlier investigators * engaged in modelling of flows 

generated by an impinging gas jet considered that the shearing forces 
acting across the cavity surfaces is primarily responsible for 
producing the bulk flow recircuJ ation. These studies, therefore, 
considered that the gas impinging over the bath surface forms a cavity 
and as the gas stream flows radially over the cavity surface, the 
outflowing gas stream drags a layer of fluid through shearing actions 
and this phenomenon produces the flow recirculation observed in such 
systems. In accordance with these concepts, therefore, the relevant 
shearing forces have been estimated from the gas injection parameters 
C viz. .section II. 2. ID and embodied in the axial and the radial 
direction momentum equations as the appropriate driving body force 
term. The relationship between the gas injection parameters and the 


to 


car-respondi ng sho’&r stress adopted in i.h«* pr©rs«*nt, study has b«N»n th& 
same as that suggested by Zhang and coworkers . The concepts 
applied to estimate the required body forces from the empirical 

c si 3 

expression of shearing stresses have been considered briefly 

alredy in section II. 4. 

To assess the effectiveness of the above concept «.e. g. .that 
the flows are essentially driven via t he -shear i ng ' forces j as applied 
to the modelling of an axi symmetric top gas injection system, an 

experimental fluid flow study reported earlier on such systems by 

c 31 y c 3i y 

Jagannathan and coworkers has been considered. These authors 

reported extensive measurements of flow parameters in a water model of 

a LD reactor via Laser Doppler Annemometer. The corresponding 

experimental conditions have been summarised in Table III. 4 for ready 

reference. The configuration of the idealised cavity surface contour 

with reference to the finite difference grid systems used for the 

present mathematical representation has been illustrated in Fig. JIT. 7. 

This shows the segmentised cavity contour , 1 he comporieni s of ‘h* 

shearing forces at any location on the cavity surface and the relevant 

normal area across which the stress components have been assumed to be 

acting. Further, the shearing forces thus estimated have been assumed 

to be prelevant over the control volumes lying immediately below the 

cavity. This implicitly implies that a shear layer thickness 

equivalent to about a control volume dimension Ce. g. , = 20mnD has been 

considered for the present flow calculations. Incorporating adequate 

expressions for F and F in relevant equations of motion In terms 

s.z s.r 

of -rCrJ and the appropriate grid parameters, and further considering 

F b =0 Csee section II . 2. 1Z> . computations were carried out corresponding 

C315 

to the experimental conditions of Jagannathan et al . 



Assumed shear layer thickness 



Shearing force along r, t(r) r rAz 
Shearing force along z,r(r) 2 rAr 


Fig. III. 7. The shearing forces acting at a location on the 
cavity surface and the configuration of the shear 
layer thickness with reference to the finite 
difference grid system. 



Numerically predicted vertical velocity components at four 
different depths in the model has been compared directly in Fig. 111. 8 
with the exper i mental observations of J aganna t han and coworkers. This 
clearly indicates considerable differences between theory and 
experiments and the predicted flow is seen to be si gni f i cantl y smaller 
C by about a factor of B or so.) than the corresponding experimental 
val ues . 

This discrepancy between theory and experiments clearly 

requires some rigorous analysis and to this end the sensitivity of 

calculations to the value of 1^ C=di stance of seperation between the 

lance nozzle and the bath free surface) must first be assessed, since 

C31 ) 

no value of 1 A has been specified by J aganna than and coworkers . It 

is important to note here that calculations reported in Fig. III. 8 has 
been driven by assigning an arbitrary but likely value of 150mm to 
1^. Similarly the assumption of a shear layer thickness equivalent to 
the size of a c on t r ol vol ume ma y a 1 so 1 nt r oduc e s ome uncer t a i ni t y in 
the numer i c a 1 pr edi c t i ons . 

To evaluate these, sensitivity of flow predictions has been 
carried out to the choice of values of 1^ and the shear layer 
thickness. In Fig. III. 9, predicted vertical velocity component at two 
different depths in the model have been compared considering two 
arbitrary but probable values of the parameter 1 Q C =1 50mm and 220irm 
respectively). This evidently indicates that predicted results are 
relatively insensitive to the precise value of 1 q chosen. Similarly, 
the shear layer thickness has been ar bitrarly increased by a factor of 
two C shear stress assumed to be pr el a vent over a region that is twice 
as thick as that considered in the computation of Fig III. 8) and the 
numerical predictions repeated for the conditions of computations 




Fig. III. 8. Numerically predicted vertical velocity component at different depths and 
their comparison with experimental measurements of Jagannathan et al. 31 
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Fig. III. 9. Sensitivity of numerical predictions to 

the values of l 0 illustrating the variation 
of vertical velocity component at two 
different depths in the model reactor . 
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cor respond! ng to Fig. III. 8. The* resul is have been ill ustr abed in 
Fig. HI. lO where the vertical velocity components at two different 
depths have been compared for two different set of values of the shear 
layer thickness. Cl early, the assumption of a shear layer thickness in 
the range Of 20-40 mm appears to be not critical to the overall 
prediction of the flow fields. These as a consequence suggest that 
uncertainities assosj ated with the values of 1 and the shear layer 
thickness are not likely to be the possible reasons for the 
discrepancy between theory and experiments illustrated in Fig. III. 8. 

111.3*2 Some additional considerations in the mathematical modelling 
of flows generated by an impinging gas jet 

The impinging gas jet in a top blown system displaces fluid 

from the cavity region and consequently, leads to an inf intismal 

increase in the bath height. Under steady state conditions when the 

bath height assumes a constant value during the injection process, it 

is readily apparent that some energy must be continuously fed into the 

system to sustain the increased bath height. Towards this, a sampie 

cal oil at i on C cor res pond i ng to t he exper i menial coridi t i ons of 

Jagannat.han and cawor ker s C ^ ^suggests that displacement of fluid from 

the cavity region necessi ales a force that is approximately equivalent 

about 13. 8N while cor r esponding contribution from the integrated shear 

—2 

force appears to be only 6.213x10 N. Evidently* the shearing forces 
alone could not have produced the increased bath height and hence 
displace an e^ui v,al ent amount of fluid from the cavity region. 
Consequently, it is legitimate to consider an appropriate gravity 
force (.mass of the fluid displaced ^ aeceJ arat i on due to gr avi tyj in 
ad»1i f j on t„c* the shearing forces, a? an additional driving force 
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Fig. III.IO. Sensitivity of numerical predictions to 
the choice of the shear layer thickness, 
illustrating the variation of vertical 
velo : it y component at two different 
depths in the model reactor. 
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producing the flow recirculation in such systems . 

On the hasi s of these concepts, at each nodal point 
immediately below the cavity surface, a representative body force has 
been evaluated C vi a . , i r.t er ms of mass of displaced fluid and the 
acceleration due t.» avi ly> added < o l he a>d al direction momentum 
equation i ri addition to the relevant shearing force components and the 
numerical predictions carried out again corresponding to the 
experimental conditions of Jagannathan and coworkers . A comparison 
between predicted and experimental vertical velocity component at four 
different depths in the model has again been illustrated in 
Fig.III.il- The incorporation of the force of impingement into the 
mathematical model is clearly seen to produce estimates of flows that 
Cin contrast to Fig. III. 8} gives excellent agreement with the reported 
experimental measurements. 

III. 4 Mathematical Modelling Of The Combination Blown Process And 
comparison Of Numerical Predictions With Experimental Measurements 

Mathematical modelling of the top as well as the bottom gas 
injection processes has been discussed extensively in the preceeding 
chapters. Since, the appropriate mechanisms of momentum transfer from 
gas to the liquid has been identified for both top and bottom gas 
injection systems, consequently, the concepts and the procedures can 
be combined to provide an effective hydrodynamic model for the 
analysis of the combination blown steelmaking process. 

To assess the ef f ecti veness of predictions via such an 
approach Cviz. , coupling the models for bottom C section III. 2. l.cD and 
top gas injection Csection III . 2. 2DsystemsD , calculations were 
performed corresponding to the exper i mewntal conditions of Zhang and 
coworkers 1 ' 215 . These experimental conditions are summarised in Table 



Numerically predicted vertical velocity component at different depths and 
their comparison with experimental measurements of Jagannathan et al . 31 


III. 4 . In Fig. III. 12 numerically predicted flow patterns in the 

combination blown model vessel has been shown. This on comparison 

with equivalent experimental measurements of Zhang et 
C 21 D 

al - Csee Fig. III. 1 3D indicate reasonable agreement between theory and 
obser vat i ons . 

In deducing the results presented in Fig. III. 12, the cavity 

volume has been artificially excluded from the domain of calculation 

by applying the ’cell porosity’ concept of Moult and coworker s C26:> . 

The principle involved has been briefly outlined in Fig. III. 14. 

There , it, is readily seen t hat. by assigning a large S val ue to "those 

P 

nodal points within the cavity no flow condition can be at once 

prescribed. Similarly, for the control volumes lying partly within the 

cavity and partly wi t hi n the bulk phase cell porosity factors Cf'sD 

have been estimated geometrically to deduce the effective convection 

and diffusion fluxes. Through this technique, the influence of the 

cavity surface and the cavity volume on the bulk phase hydrodynamics 

can be adequately taken into account. 

Towards these, it is instructive to note here that numerical 

C SI y 

procedures of Zhang and coworkers as well as that of Asai and 

C 17 } 

Szekely didnot consider the above and instead regarded the cavity 

as a part of the calculation domain. The cavity volume typically is 

C 21 D 

much smaller Cfor the experimental conditions of Zhang et al . , 

cavity volume = 0.1% of the reactor volume or cavity volume = 5% of 
the plume volume) in comparison to that of the bulk and it may 
therefore appear at the first instance that inclusion/exclusion of the 
cavity volume into/from the calculation domain may not effect the 
predicted results significantly. 








1. A completely blocked u control volume 

— ► S p :-10 30 and S u = 0.0 

2. A partly blocked u control volume 

— ► A^=A|yj , Ag = (1 - f-j)A$ , A^ = (1 Ag:, 

A W = A W J ore bl° c k Q 9 e ratios 

Fig. Ill .14. The configuration of the cavity surface with 
reference to the finite difference grid, illust- 
rating the principles involved in the cell porosity 
method applied to exclude the cavity volume from 
the domain of calculation . 


Table III. 4 Experimental conditions of different investigators 
used in the computation of results presented in Figs. III. 7 
through III. 16. 



Expert mental 

condi ti ons of 


Zhang and co— 

Jagannathan and co— 


. C213 

workers 

. C31D 

wor k er s 

Bath height , m 

O. Ill 

O. 324 

Radi as of the 
vessel , m 

0. 145 

O. 5 

Top gas flow rate. 

1 . 58x1 O' 3 

8. 33x1 0 -3 

Bottom gas flow rate. 

3. 75x1 0 -5 


3/s 

m 



I nner di ameter of 



orifice of lance 
nozzle , m 

6x10 

2. 3x10 

Distance between 
orifice of lance 
nozzle and free 
surf ace, m 

0. 154 

0. 15* 

Depth of the 
bath cavity, m 

O. 0179 

O. 08 

Radius of the 
bath cavity, m 

O. 0453 

O. 075 


Assumed 



To test computations wore- carried out again for the? 

CS1 D 

experimental conditions of Zhang et al . incorporating the cavity 

volume into the domain of calculation. In Fig. III. 15 predicted 
vertical velocity component at two different depths in the model 
vessel have been compar ed wi t h equi va i ent pr edi c t i ons i gnor i ng t he 
cavity volume. Clearly the predicted results, as one might have 
anticipated are not critical to either inclusion or exclusion of the 
cavity volume into/from the domain of calculation. It is important to 
note here that predicted vertical component of the flow deduced by 
ignoring the cavity, as seen from Fig. III. 15 is somewhat larger than 
its counterpart. This to be expected since with the exclusion of the 
cavity from the calculation scheme, the net buoyancy force acting on 
the system is accordingly also under estimated and hence this leads to 
an overall slowering of fluid motion in the system. It Is also 
reasonable to anticipate here that under such conditions if the gas 
volume fraction with in the plume is increased proportionately , the 
flows deduced from either of the calculation procedure would be 
practically equivalent. 

Finally, in Fig. III. 16, the corresponding assosiated 

turbulence kinetic energy field distribution is shown. This indicates 

-3 2 -2 

that higher turbulent regions C >10 ms 3 are restricted close to 

the free surface region as well as the region immediately below the 

cavity while elsewhere, the turbulence kinetic energy is relatively 

low. Close to the bottom as well as the central portion of the reactor 

-4 2 -2 

the turbulence is seen to be extremely low C <10 m s D . These 
together with the predicted velocity fields may provide useful 
inferences about refractory wear, scrap melting etc. in an actual 
industrial steel making systems. 



Radial distance t mxlO 

Fig. Ill .15. Theoretically predicted vertical velocity 
component at two different depths in 
the model vessel 21 illustrating the effect 
of considering /excluding the cavity 
volume on the computed results. 





A 

Fig. III. 16. Prediction of turbulence kinetic energy (m s' 1 ) 
distribution in water model of a combination 
blown reactor (corresponding to the experi- 
mental conditions of Zhang et at 21 as shown 
in Table III.4). 



Chapter IV 


CONCLUDING REMARKS 

The following major conclusions earn be drawn from the 
present investigation : 

Cl} Four different procedures to model the bottom gas 
injection phenomena Cviz. »the bulk flow, the gas volume fraction 
in the plume etc.} have been rigorously assessed and it has been 
demonstrared that the buoyancy forces imparted by the rising 
bubbles in such systems is best represented if a constant plume 
rise velocity together with bubble slip is considered in the 
mathemati cal model . 

CS} The comparisons between the theory and the experiments 
illustrate that the sheaiing forces imposed by the top gas 
injection are not likely to be the only mechanism producing the 
bath recircul ation in top gas injection systems and that a 
mechanism such as the force of impingement C estimated interms of 
the appropriate gravity force} plays a crucial role in driving the 
flow field. This concept is new and is applied for the first time 
to model an impinging gas jet induced flow system. 

C 3} It has been shown that the force of impingement is many 
folds larger than the corresponding shear forces and consequently, 
an approach based on this has been adopted to model the top gas 
injection induced motion. The predicted results thus obtained are 
shown to agree reasonably well with the experimental observations 
reported in the literature. 
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C40 Numerically predicted results also indicate that 
estimates of velocity etc. in top or combination blown system are 
likeky to be relatively insensitive to the assumed shear layer 
thickness, provided a sufficiently large number of grid points are 
applied to the numerical solution scheme. 

C53 In modelling the top gas injection induced flows, the 
cell porosity method has been applied to exclude the cavity volume 
from the calculation domain and the predicted results thus 
obtained indicate that exclusion/inclusion of cavity volume 
from/into the calculation domain have only marginal effect on the 
overall accuracy of predicted results. 

C 63 Combining the appropriate concepts for the top and the 
bottom gas injection procedures, a mathematical model for 
computation of flows in a combination blown steel making reactors 
has been developed and the numerical predictions evaluated against 
reported experimental measurements deduced from the water model 
studies and reasonable agreement between the two demonstrated. 



Chapter V 


RECOMMENDATIONS FOR FUTURE WORK 


The scope of the present work can be further expanded 
considering the following : 

Cl 3 The axi symmetric mathematical model developed in the 
present study Is to be applied to investigate mixing and 
mass “transfer phenomena in such systems. 

C B'J An idealised cylindrical vessel geometry has been 
cons j dered in the development of the present model. The model must 
be made broad enough] in its scope of computation so that an actual 
indusiriai configuration C which are slightly distorted cylinders^ 
can be employed, 
and finally, 

C3} The axi symmetric model applied in the present study can 
be modified/improved to a fully three dimensional model so as to 
study theoretically the effect of multiple tuyeres or nozzle which 
are typical of industrial combination blown systems. 
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APPENDIX 



TO FIND THE SOLUTION OF A DIFFERENTIAL EQUATION 
duCz3/dz + uCz3/z = pi *Q*g/C2C tanCl amada3 *uC z3 *z3 **23 BY 
RUNGA -KUTTA METHOD WITH THE HELP OF GIVEN BOUNDARY 
CONDI TI ONS z =zC 1 3 , u fz) =uC 03 . THE SOLUTI ON OF THE ABOVE 
DIFFERENTIAL EQUATION GIVES THE VELOCITY AS A FUNCTION 
OF z. 

OPENC UNI T=23 , FI LE= * RUNGA. IN’) 

OPENC UNI T=22 , FI LE= * RUNGA . OUT ' 3 

IMPLICIT REAL*8CA-H,0-Z3 

DI MENSI ON ZC 503 , UC 503 , ALFAC 503 , RC 503 

READC23,*3 C RC N3 , N=1 . 21 3 

ZC 1 3 =0. 01 842 

UC 1 3 =1 8 . 86 

US=0. 2259 

H=0. 003335 

NMAX =21 

PI =3. 1 4 

Q=3. 3333E-05 

G=9. 81 

TANLAM=D. 086 
DO 20 N=1 , NMAX-1 

SI =PI *Q*G/C 2*C TANLAM*UC N3 *ZC N3 3 **23 -UC N3 /ZC N3 
ZC N3 =ZC N3 +H/2 
UC N3 =UC N3 +S1 *C H/23 

S2=PI *Q*G/C 2*C TANLAM*UC N3 *ZC N3 3 **23 -UC N3 /ZC N3 
ZC N3 =ZC N3 +H/2 
UC N3 =UC N3 +S2*C H/23 

S3=PI *Q*G/C 2C TANLAM*UC N3 *ZC N3 3 **23 -UC N3 /ZC N3 
ZCN3=ZCN3 +H 
UCN3 =UCN3 +S3*H 

S4 =PI *Q*G/C 2C TANLAM*UC N3 *ZC N3 3 3 **23 -UC N3 /ZC N3 
ZC N+l 3 =ZC N3 +H 

UC N+l 3 =UC N3 +C H/63 *C SI +2*S2 +2*S3 +S43 

CONTI NUE 

DO 40 N=1 . NMAX 

A=PI *RC N3 **2*C UC N+l 3 +US3 

F=4 . 0*PI *RC N3 **2*US*Q 

P=2. 0*PI *RC N3 **2*US 

ALFAC N3 —A/P -C DSQRTC A**2 -F3 3 /P 

CONTINUE 

WRI TEC 22, 4003 

FORMATC 8X , ’ NO* , 9X , * ZC M3 * , 9X , ’ VELC m/s3 ’ , 8X , * RADC M3 ’ 

, ’ 6X , * ALFAC N3 * /3 
130 30 N=1 » NMAX 

WRI TEC 22 , 5003 N , ZC N3 , UC N3 , RC N3 , ALFAC N3 

FORMATC5X,I4,5X,F10. 4,5X,F10. 4.5X.F10. 4.3X.F10. 6/3 

CONTINUE 

STOP 

END 



